Tidal Waves in ^°^Pd: A Rotating Condensate of Multiple d bosons 
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Low-lying collective excitations in even-even vibrational and transitional nuclei may be described 
semi- classically as quadrupole running waves on the surface of the nucleus ("tidal waves"), and the 
observed vibrational-rotational behavior can be thought of as resulting from a rotating condensate 
of interacting d bosons. These concepts have been investigated by measuring lifetimes of the levels 
in the yrast band of the ^°^Pd nucleus with the Doppler Shift Attenuation Method. The extracted 
B(E2) reduced transition probabilities for the yrast band display a monotonic increase with spin, 
in agreement with the interpretation based on rotation-induced condensation of aligned d bosons. 

PACS numbers: 21.10.Tg, 27.60.+j, 21.10.Ky, 21.60.Ev 



Collective quadrupole excitations of nuclei are classi- 
fied as "rotational" and "vibrational" with the rigid rotor 
and the harmonic vibrator being the limiting ideal cases. 
The harmonic vibrations {d bosons) are five-fold degen- 
erate with respect to the angular momentum projection; 
they couple to multiplets of different angular momen- 
tum. Rotational bands that extend over ten or more 
states are ubiquitous. Vibrational excitations of spheri- 
cal nuclei, on the other hand, are much less distinct since 
the time scale of collective vibrations is not much larger 
than that of intrinsic quasiparticle excitations. The adi- 
abatic separation of the time scales, which is a prereq- 
uisite for the appearance of collective quantum states, 
rapidly deteriorates with the number of excited quanta 
(phonons). While evidence for the two-phonon triplet is 
often observed, identification of all members of the three- 
phonon multiplet is already problematic (see, for exam- 
ple, Ref. [1]). Fig. [1] schematically shows the location 
of the vibrational states and of the quasiparticle excita- 
tions for a spherical nucleus. With increasing phonon 
number n, the collective states are embedded into a pro- 
gressively dense background of quasiparticle excitations. 
The coupling to this quasiparticle background fragments 
the collective levels and they cease to exist as individual 
quantum states. The density of quasiparticle excitations 
is lowest near the yrast line, which is the sequence of 
states with the lowest excitation energy for a given an- 
gular momentum /. With increasing phonon number n, 
the yrast members of the vibrational multiplets are ex- 
pected to keep their identity as collective quantum states 
the longest. 

In this Letter, we report on the first identification of 
a seven-phonon yrast state, which can be interpreted as 
a rotating condensate of seven d bosons that align their 
spins. The yrast states of ^^^Pd are nearly equidistant 
in energy up to spin / = 16, which is characteristic of a 



vibrational sequence. Indeed, the existence of an excep- 
tionally long vibrational band in ^^^Pd was already noted 
by Regan et al. [2], making this nucleus an attractive 
case for this investigation. However, in order for these 
states to represent a sequence of stacked d bosons, the 
reduced transition probabilities between adjacent states, 
B{E2; I 2), must increase linearly with the number 

of bosons. In this work, lifetime measurements are pre- 
sented that confirm this expectation. The measurements 
are complemented by theory, which provides a connec- 
tion with the underlying microscopic structure that has 
remained a challenge in the case of multi-phonon excita- 
tions. 

Fig. [2] illustrates the vibrational modes of minimal 
and maximal angular momentum of a classical ideal liq- 




FIG. 1: (Color online) Schematic representation of the loca- 
tion of the collective quadrupole vibrational excitations rela- 
tive to the quasiparticle excitations. The darker shades ap- 
proximately indicate higher densities of quasiparticle states. 
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FIG. 2: Flow pattern of a vibrating droplet of an ideal liquid. 
The arrows indicate the motion of the surface. For a given 
energy, there are degenerate modes, which differ by their an- 
gular momentum. Left: The oscillating motion, which car- 
ries no angular momentum. Right: The traveling wave (tidal 
wave), which carries the maximal angular momentum (yrast 
state). The outward motion of the surface at south-east and 
the inward motion at north-east combine such that the east- 
ern crest moves down. 

uid. The oscillating standing wave on the left-hand side 
represents the zero-angular momentum members of the 
vibrational multiplets, which are most rapidly drowned 
in the sea of quasiparticle excitations. The vibrational 
yrast states correspond to a wave traveling over the nu- 
clear surface, as illustrated on the right hand side. The 
surface rotates with the constant angular velocity uj as 
is the case for the rotation of a rigid body. However, 
the flow pattern is irrotational. As is characteristic of 
a surface wave, liquid moves from the wave front under 
the crest to the back side. The energy and the angular 
momentum increase with the amplitude of the wave, but 
the frequency remains constant. In the case of rigid ro- 
tation, the energy and the angular momentum increase 
with the angular frequency while the shape remains un- 
changed. The name "tidal wave" has been suggested for 
the yrast mode Q in analogy to the propagation of ocean 
tidal waves. The tidal wave then corresponds to a wave 
packet of aligned d bosons, which rotates with constant 
angular velocity. The tidal wave concept has been ap- 
plied previously to near-equidistant band structures of 
different spin-parity sequences. Specifically, Refs. 0, [HI 
interpreted the A/ = 1 sequences of fixed parity in ^^^Os 
as a tidal wave running over a triaxial surface, and Ref. [g^ 
interpreted the alternating-parity sequences in 220^^]^ 
a reflection-asymmetric tidal wave traveling over a spher- 
ical core. 

In terms of a rotational interpretation, the difference 
between a rotor and a tidal wave lies in the way in which 
angular momentum is generated. A rotor generates angu- 
lar momentum by increasing the angular frequency uo at 
nearly constant deformation (and, hence, constant mo- 
ment of inertia). A tidal wave nucleus, on the other 
hand, generates angular momentum by increasing defor- 



mation (i.e., changing /3, 7) at nearly constant angular 
frequency uj. For most transitional nuclei, the situation 
is intermediate: Both mechanisms of generating angu- 
lar momentum are present, one favored over the other as 
the angular momentum increases. The rotor-like scenario 
has been well studied, with the measured reduced transi- 
tion probabilities, i.e., B{E2) values, confirming that the 
deformation remains rather constant [7]. However, the 
strong increase of deformation with angular momentum, 
which characterizes a tidal wave, is demonstrated for the 
first time in the work reported here. 

The mean lifetimes of yrast states in ^^^Pd were 
measured using the Doppler Shift Attenuation Method 
(DSAM). The experiment was performed at the Argonne 
Tandem Linear Accelerator System (ATLAS) facility us- 
ing the ^^Ge(^^Si, 4n) reaction, at a beam energy of 110 
MeV. The target consisted of 500 /ig/cm^ isotopically 
enriched ^^Ge backed by a 26 mg/cm^-thick layer of Au; 
the backing thickness was sufficient to stop the recoiling 
nuclei. A thin layer of 11 /ig/cm^ Al was sandwiched 
between the target and backing to prevent the migration 
of Ge into Au. The emitted 7 rays were detected with 
the Gammasphere array ^] which, at the time of the 
experiment, comprised ninety eight Compton-suppressed 
high-purity Ge detectors, arranged in 16 rings of con- 
stant angles relative to the beam direction; a total of 
3.1x10^ three- and higher-fold coincidence events were 
accumulated. The low-spin structure of ^^^Pd is already 
well established (9l-[ll|. The present investigation has fo- 
cused primarily on the measurement of the lifetimes of 
members of the yrast band. 

The accumulated data were sorted angle-by-angle into 
several unique- fold event databases, using the BLUE pro- 
gram [12|. For each database, an angle-dependent back- 



ground subtraction algorithm was applied, and the 
resulting, background-subtracted, databases were double 
gated in order to extract coincidence and error spectra 
for each ring of Gammasphere for further analysis. 

Doppler-shifted and broadened lineshapes were ob- 
served for transitions in the yrast band up the 16+, 
7244-keV level and lifetime analyses were performed us- 
ing a modified version of the LINESHAPE analysis pack- 
age [ij. A total of 50000 Monte Carlo simulations 
of the velocity histories of the recoiling nuclei travers- 
ing the target and backing materials were generated in 
time steps of 0.0089 ps and were subsequently converted 
into time-dependent velocity profiles. The total number 
of time steps was set to 198 and the time step value 
was determined in accordance with the discussion in 
Ref. [14]. Electronic stopping powers were taken from 
Ziegler's tabulation [l5|, with low-energy modifications. 
Energies of in-band transitions were determined from fits 
to the experimental data and the side-feeding intensi- 
ties obtained directly from the measured 7-ray intensi- 
ties within the yrast band. The side feeding into each 
level and to the top of the yrast band was modeled by a 
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TABLE I: Measured lifetimes and electromagnetic transition 
probabilities for the yrast band in ^°^Pd. Experimental un- 
certainties were derived from the behavior of in the vicinity 
of the best-fit parameter values. 



E^(MeV) 


If ^ 




t(ps) 


B{E2){e%^) 


Ot(eb) 


0.836 


6+^ 


► 4+ 


1.149(32) 


0.174(22) 


2.36(15) 


0.902 


8+^ 


6+ 


0.720(42) 


0.190(11) 


2.41(07) 


0.980 


10+- 




0.409(28) 


0.221(15) 


2.56(09) 


1.062 


12+- 


10+ 


0.252(18) 


0.240(14) 


2.65(08) 


1.083 


14+- 


-> 12+ 


0.201(11) 


0.272(15) 


2.80(08) 


1.104 


16+- 


-> 14+ 


0.180(21)^ 


0.277(14) 


2.81(07) 



*The lifetime of the 16+ level could not be separated from the 
side-feeding lifetime. The value given is, therefore, the lower limit. 



five-state rotational cascade with independently variable 
lifetimes. A minimization for the observed Doppler 
shifted lineshapes was then performed using the in-band 
and side-feeding lifetimes, background, and contaminant 
peak(s) as input parameters. Experimental uncertain- 
ties in the extracted lifetimes were determined based on 
the behavior of the the vicinity of the min- 

imum t^, by a statistical method using the MI- 
NOS [l8| program. Although side-feeding effects were 
initially included in the analysis, the extracted lifetimes 
were found to be insensitive to changes in the side-feeding 
intensities. More details on the fitting procedure can be 
found in Refs. [16, 19]. 

Experimental data and lineshape fits at both forward 
and backward angles for two representative transitions 
(1104 and 1062 keV) in the yrast band are displayed in 
Fig. [3l The inserts are expanded views of the corre- 
sponding lineshapes. The extracted lifetimes of states 
in the yrast band, the corresponding reduced transition 
probabilities, and the associated transition quadrupole 
moments Qt are summarized in Table HI The reduced 
transition probabilities were obtained from the measured 
lifetimes using the expression: 



B{E2) = 



0.0816 



E^{E2)r{E2) 



(1) 



with the 7-ray energies, E^^ in MeV, and the partial life- 
times of the transitions, r{E2)^ in ps. 

Fig. |4] provides the experimental moments of inertia 
J = I/uj = 2I/[E{I) - E{I - 2)] as a function of spin 
/. The harmonic limit for the free bosons is displayed by 
the dashed line marked FB and corresponds to a constant 
LjJ equal to one half of the vibrational frequency Q. The 
experimental J7 moment is a nearly linear function of / as 
indicated by the dashed line marked IB (for interacting 
bosons). It deviates from FB by the small offset at / = 0, 
which is a measure of the anharmonicity. The angular 
momentum increases due to the increase of J7 while u 
remains nearly constant. 

The measured B{E2) values (Fig. [5]) increase lin- 
early with /, such that the ratio B{E2)/J is constant 



within the experimental uncertainties (see the insert). 
This demonstrates that the angular momentum gain orig- 
inates from the increase of the wave amplitude (defor- 
mation) while the rotational frequency does not increase 
significantly. This is in stark contrast to a typical rotor 
where the transition probabilities remain relatively con- 
stant over a given spin range, while the increase of the 
rotational frequency translates into a gain of angular mo- 
mentum. This provides the first evidence that the yrast 
line of ^^^Pd corresponds to a slightly anharmonic tidal 
wave and the measurements reported here establish the 
characteristic increase of the tidal wave amplitude up to 
/ = 14. 

It is remarkable that the anharmonicity shows up in a 
simple way as a constant upshift of the functions J {I) 
and B{E2; I 1—2). This may be interpreted as follows: 
The yrast line of ^^^Pd is characterized by an anharmonic 
wave that starts out with a small deformation (i.e., non- 
aligned d bosons are present with a small probability) 
and increases in amplitude along the sequence of states 
by adding aligned d bosons. 

The transition quadrupole moment measures the am- 
plitude of the tidal wave. For 46^Pd56, it is = 2.8 eb at 
/ = 14, to be compared with Qt = 2.96 eb for the 1 = 2 
state in the nucleus 46^Pd64, which is located in the mid- 
dle of the neutron shell [21]. In the case of ^^^Cd, the ex- 
perimental amplitude of the tidal wave reaches Qt = 2.53 
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FIG. 3: (Color online) Lineshape fits to the 1104- and 1062- 
keV transitions in the yrast band in ^°^Pd. The upper and 
lower panels correspond to the backward and forward detec- 
tors, respectively. The inserts are expanded views of the line- 
shapes. 
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FIG. 4: (Color online) Moments of inertia, as a function 
of spin, /. The dashed line marked FB is the limit for har- 
monic bosons. The dashed line marked IB indicates the near- 
linear trend of interacting bosons. WS-SCTAC represents the 
cranking + micro-macro calculations described in the text. 



eb for / = 8 [21^. [All quoted Qt values have been ob- 
tained from the lifetime information provided in the cited 
references.] In this case, the s configuration with the 
rotational alignedpair of /111/2 quasineutrons becomes 
yrast for / > 12 [2l]. Its nearly constant smaller defor- 
mation {Qt ~ 1.7 eb) thereafter signals the transition to 
the rotational regime (c.f. [3]). In the case of ^^^Pd, the 
rotational aligned s configuration does not become yrast 
until after 7=16, which, as mentioned previously, is one 
reason why this nucleus was deemed very attractive from 
the point of view of studying the boson condensate. Un- 
fortunately, it was not possible to determine the lifetimes 
of these "rotational" states from the present data. 

The tidal wave concept allows for a semiclassical calcu- 
lation of the energies and B{E2) values of the yrast states 
in vibrational and transitional nuclei. The tidal wave has 
a static deformed shape in the co-rotating frame of ref- 
erence. This has led to the microscopic description sug- 
gested in Ref. [3j , which is based on the cranking model 
without resorting to the small amplitude approximation. 
Ref. applies the shell-correction tilted-axis cranking 
(SCTAC) version [22'] of the cranking model to even-even 
nuclides with 44 < Z < 48 and 56 < < 66. In this 
approach, the energy surface is calculated by the micro- 
macro method, subject to the constraint that the expec- 
tation value of the angular momentum operator equals /. 
The energy is minimized with respect to the deformation 
parameters P and 7. Deformed solutions are found in Ref. 
[3j for / > 2, even when the solution was spherical for 
/ = 0. These calculations reproduce the collective yrast 
states rather well. They also describe the intrusion of the 
aligned /111/2 two-quasineutron states into the yrast line, 



which causes the backbending phenomenon seen in most 
of the nuclei studied. The details of the present calcu- 
lations are the same as those of Ref. [3|, with the ex- 
ception that the modified oscillator potential is replaced 
by a deformed Woods-Saxon potential with the so-called 
"universal" parameters [23]. The results are labeled by 
SCTAC in Figs, g] and [5] ( "TAC" refers only to the code; 
the self-consistent solution actually rotates about the in- 
termediate principal axis of the slightly triaxial shape). 
It should be noted that there are no parameters adjusted 
to the experiment. 

As seen in Fig. [H the observed behavior of J7 as a 
function of spin / is reproduced well by these calcula- 
tions. The calculated J7(/) moment remains close to the 
IB line because the ground state configuration has been 
employed throughout. The calculated B{E2] 1^1 — 2) 
values (Fig. [5]) also follow the experimental data overall, 
albeit with somewhat larger fluctuations than those ex- 
hibited by the experimental values. The reason for this 
apparent discrepancy lies in the fact that the deforma- 
tion parameters are those associated with the minimum 
of the energy calculated by the micro-macro method. 
This method neglects zero point fluctuations of the shape 
which, when properly accounted for, are expected to re- 
sult in washing out of these fluctuations. Thus, the semi- 
classical cranking+micro-macro calculations describe the 
yrast states of ^^^Pd up to the seven boson state rather 
well, once the anharmonicities are correctly accounted 
for. 



35 
30 



I — I — ' — I — ' — I — ' — r 
• Exp 

•A" WS-SCTAC 



—I — T" 
IB y'' 




FIG. 5: (Color online) Reduced electromagnetic transition 
probabilities for the yrast states of ^°^Pd as a function of 
the s pin , /. Lifetime data for the 2^ and 4^ states are from 
Ref. ^]. The dashed line marked FB is the limit for har- 
monic bosons. The dashed line marked IB indicates the near- 
linear trend of interacting bosons. WS-SCTAC represents the 
cranking + micro-macro calculations described in the text. 
The insert provides the experimental and calculated ratios 
B{E2)/J as a function of /. 
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In summary, by measuring the lifetimes of the yrast 
states in ^^^Pd up to a spin of / = 14, the first clear 
evidence has been provided for rotational-induced con- 
densation of aligned d bosons. The linear increase of 
the reduced transition probability with the boson num- 
ber n = //2, which is expected for condensation, was 
observed up to n=7. The mutual interaction of bosons 
causes anharmonicities which appear as a constant up- 
shift of the reduced transition probability and the mo- 
ment of inertia as function of spin /. This upshift may be 
interpreted as caused by a small fraction of non-aligned 
d bosons to which the aligned d bosons are added. Semi- 
classically, the condensate represents a tidal wave trav- 
eling over the nuclear surface with constant angular ve- 
locity (equal to one half of the vibrational frequency), 
where the angular momentum gain arises from the in- 
crease of the wave amplitude. Since the wave motion is 
not quite harmonic, a slight increase of the rotational fre- 
quency with spin is seen. The tidal wave is described in 
the framework of the cranking model based on the micro- 
macro method, which describes the data well without any 
adjustment of the model parameters. 
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